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ABSTRACT 


A theoretical investigation of the raicroatructure of solar wind stream 
interaction regions is presented. We discuss the role of several electro- 
static kinetic instabilities which may be important within the stream 
interface and the compression region. Inside of 1 AU the interface is 

likely to be stable against the electrostatic streaming instabilities 
considered. Between 1 and 2 AU we argue that the interface will excite thf 
magnetized ion-ion instability. The compression region is also found to be 
unstable beyond 1 AU where the modified two-stream instability, beam- 
cyclotron instability, and ion-acoustic instability will be important in 
determining the structure of the compressive pulses as they evolve into 
forward and reverse shocks. We conclude that the modified two-stream 
instability and beam-cyclotron instability predominately play a role in 
heating the electrons to the threshold for the ion-acoustic instability. 
Various electrostatic plasma waves, ranging in frequency from the lower- 
hybrid to harmonics of the electron cyclotron frequency, would be produced 
by these instabilities . Their signature should also be seen by high time 
resolution measurements of the temperature of the various plasma species. 
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1. INTRODUCTION 


Two prominant features of the interaction regions of high and low speed 
solar wind streams are the stream interface and the surrounding compression 
region which, at distances beyond 1 AU, often evolves into forward and 
reverse shock pairs. In both these regions there is observational and 
theoretical evidence suggesting that strong gradients form _n solar wind 
velocity, plasma density and/or magnetic field. In this paper, we examine 
the consequences of a.ssuming that as the scale lengths of the gradients 
decrease both the stream interface and compression region will excite 
several electrostatic instabilities. 

In the following section we first investigate the stream interface, and 
argue that ions from the fast stream are likely to penetrate at least an ion 
Larmor radius into the slow stream, thereby exciting the magnetized ion-ion 
instability. The region in the heliosphere where this instability can be 
excited is restricted by the necessity for there to be a fairly large (>20 
km s”^) discontinuous jump in the component of solar wind velocity ortho- 
gonal to the interplanetary magnetic field. This is most likely to first 
occur near 1 AU or slightly beyond. The wave modes excited, their fre- 
quency, and the effect of the instability on the temperatures of the 
electrons, protons and helium are discussed. 

In §3 a similar discussion is developed for the compression region. 
There too, we argue that the plasma is likely to become unstable as the 
forward and reverse shocks begin to forvi. The electron-ion instabilities 
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excited there will produce a wide range of electrostatic waves, ranging in 
frequency fs'om the lower-hybrid to harmonics of the electron gyrcfrequency . 
Estimates are given for the saturation amplitudes of these instabilities, 
and their role in limiting the thickness of the developing shook structure. 
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2 . STREAM INTERFACES 


Observations. A stream interface is a sharp transition in velocity, 
density, temperature and flow angle of the wind often observed within the 
stream interaction region. It separates the cool dense slow stream from the 
hot tenuous fast stream and is associated with a local maximum in thermal 
pressure [for a review see Burlaga, 1975]. The existence of a tangential 
discontinuity separating the fast and slow streams was suggested as early as 
1963 [v., e.g. Dossier and Fejer, 1963]. 

Interfaces were first reported by Belcher and Davis [1971] based on 
three-hour average data. Subsequently Burlaga [1974] showed that interfaces 
could be very thin although he also included fairly broad interfaces where 
the transitions in flow velocity and density took -30 minutes, but these 
will not concern us. The existence of large jumps in flow velocity at 
stream interfaces appears confined to observations near 1 AU. Beyond 2 AU 
Pioneers 10 and 11 did not usually observe such discontinuities [see, e.g. 
Smith and Wolfe, 1977]. The kinetic effects discussed below can only be 
important at those interfaces which exhibit sharp, discontinuous changes in 
flow velocity over the fastest time scales available to plasma instruments. 
One should keep in mind, however, that the time interval over which those 
measurements have traditionally been made is still long compared with the 
few tenths of a second that it takes a region of a few ion Larmor radii to 
be convected past a spacecraft. Additional observations have been reported 
by Gosling et al . [1978] who analyzed the properties of the flow near 23 
discontinuous interfaces with thicknesses less than 4x10" km. These were 
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observed only when the solar wind speed was less than 450 km . They 
further found evidence of a sharp shear flow in solar wind velocity at 
interfaces, and noted that at interfaces the electron temperature rose 
sharply there by about 40it. We will argue below that several of these 
features are a natural consequence of kinetic instabilities operating in 
regions of sharp gradients in plasma velocity ana density. 

Because of the very high conductivity of the solar wind it is difficult 
for material on neighboring magnetic flux tubes to interpenetrate. Conse- 
quently, as fast plasma overtakes slow plasma, a compression region forms 
which tends to steepen with increasing heliocentric distance, eventually 
forming a tangential discontinuity separating the compressed ambient plasma 
from the fast stream [Hund hausen, 1972, p. 132ff; Hundhausen and Burlaga, 
1975 J. If this boundary becomes as thin as a few proton Larmor radii, the 
faster stream can begin to penetrate the slower one. It .is at this stage 
that kinetic effects become important. Inside of •>-0.3 AU, where adjoining 
flux tubes are nearly radial, the only kinetic interaction allowed is a 
viscous shear [e.g. Eviatar and Wolf, 1968]. By the time the streams have 
reached 0.5 AU, however, the component of the velocity jump across the 
interface perpendicular to £ can exceed t.ie thermal velocity of solar wind 
protons. This situation can then be unstable bo excitation of an electro- 
static instability known as the magnetized ion-ion instability (Mil), pro- 
vided that 2V, < U*sin(r) < 2. 5(1+6 )V where U is the magnitude of the 

X G A 

velocity change at the interface, is the (garden-hose) angle at 

heliocentric distance R between the radius vector and the magnetic line of 

force, n =2.7x10 ^ is the angular rotation rate of the sun, 6 ^SrNT /B' is 
s e e 

the ratio of electron thermal energy density to magnetic energy density, T^ 
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•1 

is the electron temperature, V^£B(4nNM)~ ' " is the Alfven speed, N is the 
total density in the interpenetration region of the two streams, is the 
solar wind velocity in the slow stream, M is the ion mass, and is the ion 
thermal velocity [Papadopoulos et al., 1971]. The geometry of the interac- 
tion is shown in Figure 1. 

The.Ma gnetlzed Ion-Ion Instability. This instability has been investi- 
gated in a somewhat different context with regard to the stream interaction 
in papers by Papadopoulos [1973a] and Papadopoulos et al. [1974]. Our 
application of the magnetized ion-ion instability to the stream interaction 
problem differs in several respects. Papadopoulos [1973a], following 
Papadopoulos et al. [1971], simplified the dispersion relation for the Mil 
by assuming that the ions were cold, 1 .e . << w/k ~ U. In addition, as 
noted above, U must be less than a few times V^. In the solar wind is 
often of order and these two conditions are difficult to satisfy simul- 
taneously. This has necessitated a return to the more general dispersion 
relation in which the ions are warm. From the studies of Burlaga [1974] and 
Gosling et al. [1978], we know that discontinuous jumps in the solar wind 

velocity observed at 1 AU are typically only 15-50 km s~^ } and even with a 
_1 

jump of 50 km s the component directed across B_ is reduced by a factor of 
sinr . In this situation the plasma is very close to marginal stability and 
the warm ion dispersion relation must be used. 

Furthermore, Papadopoulos [1973a] and Papadopoulos et al . [1974] 
assumed that even close to the sun the plasma flowed across [cf. Figure la 
in Papadopoulos, 1973a]. But to date there is no observational evidence in 
support of such an initial condition. It is interesting to note, however, 
that the interaction regions of streams do appear sharper at smaller helio- 
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centric distances [Rosenbauer et al., 1977J, although the scale of the 
velocity gradients remains larger than required for interpenetration. The 
assumption that gradients associated with stream interaotions should be 
steeper near the sun than at 1 AU was used in the multifluid simulations 
reported by Papadopouios et al . [197^], and was implicit in the kinetic 
calculations of Goldstein and Eviatar [1973] and Papadopouios [1973a]. 
Recently, evidence has been presented indicating that the interface itself 
is not necessarily a discontinuity inside ~0.6 AU [Schwenn, Muhlhauser, and 
Marsch , 1978]. In the following discussion we adopt the picture illustrated 
in Figure 1 in which the fast stream originates in regions of the corona 
adjacent to those giving rise to slow flow; as would occur at coronal holes. 
Thus close to the corona the flow is initially characterized by a shear in 
velocity, and only at distances greater than 0.6 AU will the component of 
velocity across £. exceed 2V^. 

In deriving the dispersion relation for the Mil we have retained terms 
which arise from electromagnetic and finite effects because in the solar 
wind those terras can be important [Wagner et al., 1971]. 

The properties of the magnetized ion-ion instability can be found from 

the general dispersion relation for waves propagating at a large angle to £ 

[Stix, 1962]. In the appropriate regime of unraagnetized ions but magnetized 

electrons, the general dispersion relation simplifies greatly. These 

simplifications are justified so long as the growth rate, y, exceeds 

n^=eB/Mo, the thermal ion Larmor frequency (M is the proton mass), which in 

-1 - 1 

turn z-’equires that the wavenumber, k, satisfy << k << , where 

1/L. = P. /V. , and where L. (L ) is the thermal ion (electron) Larmor 
radius, and n^=|eB/mc| (m is the electron mass). The sums over Bessel 
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functions appearing in the general dispersion relation can then be elimina- 
ted. The frequency range is thus restricted to fl* << |(uj* « tJ*, where to is 

X 6 

the complex frequency of the oscillation. We retain the warm electron term 

in order to investigate the validity of the usual assumption that the waves 

propagate so nearly perpendicular to ^ that o)/{kV sine)»1, where tane = 

e 

k||/k^. 

With these assumptions the dispersion relation for the magnetized ion- 
ion instability in the frame of the slow plasma becomes (Appendix A) 


atij; 


ak'V* 




'u-kUsin( r)cos( 0 ) 
~/ 2 kv' 


) 


(l-Tl)!!)? ZY_ — \ ■ 

\/2kV / 

2k 


/ 

\/2kVgSin(e)/ 


where w? = ilttNe */(M,m) , T. = mV*, and a is the fractional density in each 

X I 6 0 6 

stream. As pointed out by McBride and Ott [1972], in reference to the 

related modified two-stream instability, electromagnetic effects (w*/k=‘c*^0) 

0 

can sometimes have a stabilizing influence. In the solar wind, because 
L./L ~20 one has uiVk^c^il and the electromagnetic corrections cannot be 

X 6 0 

ignored. Observations indicate [Burlaga, 197^ and Gosling et al., 1978] 
that at stream interfaces the faster plasma is about half as dense as the 
slow material, and so we nave taken a=1/3. Equation (2.1) can then be 
solved numerically. The existence of unstable roots depends sensitively on 


1+ UJ* 

e 

k*c*(1+B ) 
6 


1 + to* 

e 

n* 

e 


1 + 


k*c*( 1+6_) 
e 
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the variation of solar wind parameters throughout the interplanetary raedlUiii. 

Close to the sun, where U’sinr/V^ << 1, no instability is found. As one 
moves out in heliocentric distance to about 0.8 AU, U»slnr/V^ s 1 an*’ 
excitation of the magnetized ion-ion instability depend.*" primarily on the 
local values of U, V^, and We have investigated solutions to (2.1) 

extending from 0.6 AU to 5 AU using a range of values for solar wind 
parameters. Because the Mil excites waves propagating perpendicular to 
it is (Tj)j^ which controls the onset of the instability, and so we have used 
= 2 eV (2x10'’ °K) in (2.0 at 1 AU. For N and £ (also at 1 AU) 
we took N=7cm ^ and £=5x10“^ G. To evaluate (2.1) at heliocentric distances 
Doth larger and smaller than 1 AU, we adopted the following scaling laws: 
N(R)-1/R*, T.(R) - R ^^^3)^ .p ^ P (1/3) [3ee^ e.g. the discussion in 
Hundhausen , 1972 and the ro :ent ra • ;.'?urements by Ogilvie and Scudder, 1970]. 
To determine the variation in the threshold conditions of the magnetized 
ion-ion instability we varied from 250 km s“^ to ^50 km s~^ using the 
T-V relation [Burlaga and Ogilvie, 1973], This was normalized so 

that (Tp)j^ = 2 eV for = 350 km s"\ 

As an example, consider the solution to (2.1) at 1.5 AU with V = 300 

O n 

_ 1 

kjD s . Using the 1 AU values and the scaling relations defined above, we 

found unstable rocts to (2.1) when the discontinuous jump in velocity at the 

interface was U = 1<0 km (Recall that U denotes the total jump in 

at the interface — the component perpendicular to B is U*sinr(fl)]. In Figure 

2, we have plotted the real and imaginary roots, u)(k,6) and Y(k,0) first 

against 0 using the value of k which maximizes > at each 0 (denoted k ) and 

m 

then against k using the value of 0 which maximizes y at each k (denoted 
e^). y is plotted for two different values of Figure 2a illustrates 
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the well known property of tne Mil that it is most efficiently excited 
across the magnetic field. JJote that for 8 < 1®, we could have assumed that 
the electron.'? were "cold” and expanded the electron plasma dispersion 
function rPapadopoulos , 1973a]. The effect of making such a cold electron 
approximation Is illustrated in Figure 3 (using somewhat different para- 
meters) . The cold electron approximation is seen to be valid only so long 
as e < 1”. 

From Figure 2, we can define v =Y(k ,6 ) as the maximum value of y, 

m m m ’ 

maximized with respect to both k and 0. In Figure 4, we have plotted and 

tOjjj [the value of (u(k,e) associated with y^] against R using two different 

values of U and a wind velocit'/ of = 300 km s~\ Also plotted in the 

Figure is n. and u defined by 
1 o 


b) 


2 

O 


[ 1 + (o)V(}^) ] 
e e 


( 2 . 2 ) 


In the solar wind because wyn*»1, u =ui =/n p . Papadopoulos et al. [1971] 

c s O L»n G i 

found that Yj^ was as large as However, because of the stabilizing 

effects of the electromagnetic contributions to (2.1), together with the 
fact that U is not much greater than V^, we find that Yj^ is systematically 
smaller than u>^, but that “ “o’ cold ion approimation . 

(Note that our definition of differs by a factor of /2 from that of 
Papadopoulos et al. [1971].) For U = 50 km s”^ the magnetized ion-ion 

instability is first excited at 1 AU, while for U = 40 km s~^ threshold for 
the Mil is not reached until 1.5 AU. This reflects the distance at which 
U'sinr first exceeds by an amount sufficient to excite the instability. 
Our assumption that the ions are demagnetized is verified a posteriori 
because y^ > 
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The heliocentric distance at which the Mil is likely to reach threshold 

is a sensitive function of U and V . This relationship Is shown in Figure 

o w 

5, where U is plotted against R f'or various values of V The curves 

o w 

represent the values of R and U for which Yjjj exceeds and suggest the 

heliocentric distance at which the Mil is likely to become important. For 

example, for discontinuous velocity Jumps of 40 km s in a 250 km s solar 

wind (recall that refers to the slow stream), the Mil can be excited at 

1 AU and beyond. However, if the "slow" wind has a velocity of 450 km s"\ 

then with U=40 km s”^ the Mil is not likely to be excited inside 2 AU. In 

general, the smaller the velocity jump and the higher the wind speed, the 

further out one must go before the conditions for excitation of the Mil can 

_ 1 

be satisfied. However, even U « 20 km s is adequate to excite the Mil 

near 2 AU with ^ 250-300 km s~\ Conversely, velocity Jumps in exces i 
of 40 km s ^ are necessary in order to excite the Mil Inside 1 AU. Thus 
this instabilj is likely to be important where is low and U - 20-40 km 
s-\ 

■Saturation Effects.. The nonlinear stages of the magnetized ion-ion 

instability have been intensively investigated both theoretically and hy 

means of computer simulations [Papadopoulos et al., 1971]. The waves first 

evolve as described by quasllinear theory until ion trapping becomes 

1 /P 

important. Stabilization then occurs after a time = (M/eE*k) =• 1 /y 
(E* is the rras value of the electric field of the growing wave [Manheiraer , 
1971]). Subsequently, the ions are reflected in the potential wells of the 
waves and the waves will stop growing after being further amplified by 
another factor of e = 2.7 [Papadopoulos et al . , 1971]. From the definition 
of T,. the maximum value of the electric field is 
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(2.3) 


= [87r(2.7)Mc^N]^^^ (k„c/u.) {y/kc)^ 

M m 1 'm m 

In Figure 6 we have plotted Ej^ between 1 and 5 AU for Vgy=300 km s" ' 
and U=50 km s~^ . Values lie are between 100-^00 uV/m with the peak occur- 
ring near 2 AU. Observation of these waves has thus far not been reported, 
perhaps because the frequency range is below that of most experiments that 
have been flown. For example, on UsHos 1 and 2 where u, „/2n < 10 Hz the 
plasma wave experiment was insensitive below 31 Hz [Gurnett and Anderson, 
1977). The waves may be difficult to detect because the instability occurs 
over only a few ion Larmor radii - 50 km at 1 AU), a distance which is 
convected past a spacecraft in less than a second. 

Stabilization by ion trapping suggests that within the interface the 
ion temperature should be significantly enhanced over its value on either 
side. In essence, within the stabilization region the thermal velocity will 
have increased to 'll [ V^+U*sin( r )] , so that as long as U’sinr 2V^^ the 
increaue in temperature can exceed a factor of two. Computer simulations in 
which U*sinr>>V^ have found that the ion temperature at stabilization is 
- 0.3Wj^[U*sinr] ^ [Papadopoulos , 1973a] i with the increase being primarily in 
(T^)j^. Gosling et al. [1978] using one-hour averages reported an increase 
of nearly a factor of two in (T^)j^ at the interface (cf. their Figure 1<) . 
With five-minute averages only the change in (T^) j | was reported [v. Figure 
2 of Gosling et al., 1978]. Of course, enhancements in caused by kinetic 
interactions are in addition to those expected due to the compression of the 
solar wind fluid at the interface and those resulting from the different 
boundary conditions in the coronal source regions of the fast and slow 
streams . 
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Alpha particles will also be affected by the Mil. Once therraalized, 
the alpha particle temperature will be of order T - (M /2) * (U*altir)*, so 

a u 

that = 4. Gosling et al. [1978] only report one-hour 

averages of alpha particle parameters near interfaces, and 30 their observa- 
tions may primarily reflect fluid, as opposed to kinetio, oharacteristios of 
the stream interaction. Nonetheless, they found that at interfaces the 
alpha to proton temperature ratio increased by more that 30J, and that T 

a 

itself rose by nearly a factor of three. 

One additional consequence of this instability is that at the interface 
the amount of shear at the tangential discontinuity will be reduced because 
the primary effect of the ion trapping is to thermalize that component of 
the flow directed across i.e. the fraction (1/2)M(U*sinr)^ . It is impor- 
tant to realize that the instability does not convert all of the streaming 
energy into electrostatic waves. In fact, for the values of shown in 
Figure 6, only a very small fraction of the streaming energy (< is 

converted into waves in the process of therraallzing the ions and producing 
the increase in proton temperature at the interface. Because the component 
of the flow parallel to fl. is unaffected, the directed flow of the fast 
plasma appears to change from radial to shear at an Interface. At greater 
heliocentric distances & larger fraction of the flow energy will be thermal- 
IzeJ by the InstablJ.ity , will increase, and the shear should decrease 
until the instability has ercied the velocity gradients to the point where 
it cannot be excited. In this context it is not. surprising that discon- 
tinuous jumps in velocity are often absent in interfaces observed beyond 1 
AU. 

ELectron Heating. At first glance the magnetized ion-ion instability 
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would not appear to be a very efficient means of heating electrons. How- 
ever, the electrons will respond to the electrostatic fields produced by the 
ion-ion interaction by drifting in the direction of £x£ [Papadopoulos et 
al., 1971]. For electi’ic fields approaching 1 mV/m, the drift velocity can 
exceed 100 km s \ which is adequate, as we shall discuss in more det' 1 in 
the next section, to excite the modified two-stream instabil y. his 

instability both heats the electrons, increasing (T ) ■ i , and lowers U , 

e 1 1 e 

their mean drift velocity. The primary source of free energy is the 
electron drift energy (1/2)raU^ [Lampe et al. , 1975]. Assuming that U = 

C 0 

- 1 - 1 

200km s , and « 1000 km s , the increase in perpendicular electron 
temperature at the interface will be - 5 %. Gosling et al. [1978] utilizing 
five-minute averaged data do report an increase in (T ) . , of some 10% at 
interfaces (v. their Figure 3). 

In the following section we examine several kinetic interactions which 
may be important in the formation of the forward and reverse shocks often 
observed in stream interaction regions beyond 1 AU. 
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3. FORWARD AND REVERSE SHOCKS 


The compression region on either side of the stream interface steepens 
with heliocentric distance until a forward and reverse shock pair forms (see 
Fig. 1). These corotating shock pa^rs have not been observed insidj 1 AU 
[Rosenbauer et al., 1977], and are only rarely observed at 1 AU [Ogilvie, 
1972]. However, when Pioneer 10 and 11 encountered corotating stream 
interaction regions beyond 1 AU, observation of shock pairs became common 
[Smith and Wolfe, 1977]. Plzzo [1978a] and [1978b], using multidimensional 
fluid models of corotating streams, has shown that formation of forward and 
reverse shock pairs is expected to occur only near and beyond 1 AU. 
Nonradial flows in these models transport mass, energy and momentum away 
from the compression region and that tends to delay the formation of shocks 
to beyond 1 AU. 

In this section we discuss several microinstabilities likely to be 
excited as the forward and reverse shocks form in the compression region. 
The instabilities considered are not meant to exhaust all possioilities , and 
for example we consider neither the closely related lower-hybrid drift 
instability [Davidson et al., 1977], nor the drift-cyclotron instability 
[Gladd and Huba, 1979]. The ones treated do, however, illustrate the range 
of frequencies and classes of waves that one might observe with appropriate 
instrumentation . 

Modified Two-Stream Instability (MTSI) . As the compressive MHD waves 
steepen, the local magnetic gradients within the waves produce streaming of 
electrons through ions. This situation was first described by Fredricks 
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[ 1969]. In a plasma in which T < 10T» this strec^ining is known to excite 

0 * 

the modified two-stream instability [Ott et al., 1972 and McBride et al., 

1972]. This instability is similar to the magnetized ion-ion instability 

discussed above; the major distinction is that instead of having two ion 

distributions flowing through each other in the presence of electrons, one 

now has a single ion distribution streaming through electrons. The MTSI can 

be excited so long as the relative velocity of electrons and ions, 

exceeds both and the sound speed, C^=(T^/M) . Again we retain the 

electromagnetic corrections to the dispersion relation, allow both the ion 

and electrons distributions to be warm, and restrict our attention to the 

-1 -1 

wavenumber and frequency range <<k<<L^ and fi|<<( oj | respectively. 

The resuiuliag dispersion relation becomes (Appendix A) 


/W-kV .G0S( 0 )\ /hi \ 

i — ^ A 

\ /2kV^ / \/2kVgSin(6)/ 


2k*V| 


2k"V^ 

e r . 


1 + 


k*o"(1+6„) 


1 + w* 
e 


e 


1 + 


k*c^(1>6 ) 
e 


(3.1) 


Excitation of the instability is determined by the magnitude of V^, which in 
turn is determined by the drift current, J , within the MHD wave. Thus 
[Fredricks, 1969] 


c(AB/Ax) = 47 tJ = 4irNeV. 

e d 




We solve (3.1) using as a free parameter. Once having found the 
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value of which maximizes the growth rate, (2.2) can be used to estimate 

the scale size of the MHD wave required to excite the MTSl . Because the 

growth rates are relatively insensitive to and R, we first solve (3.1) 
using parameters typical of 1 AU; a choice which permits comparison with the 
results of Lemons and Gary [l977Ji who have used a more general formulation 
of the MTSI to investigate the formation and structure of the earth's bow 

shock. 

In Figure 7 we plot the real and imaginary roots of (3.1) first versus 

k [at the value of e which maximizes then versus e [at the value 

of k which maximizes y(o)]. The computations were performed using = J(50 

km s"\ R = 1 AU, = (T^)j^ = 3-3 eV, B = 5 y. N = 7 cm"^, and = 0.6 and 

1.1. In this example, = 5V^. As illustrated in the Figure, as B^ 

decreases, increases as more free energy becomes availai'i. to drive the 

instability, i.e. increases. In deriving (3-1) it was necessary to 

assume that k,,/k < V./V in order to reduce to a single term the summations 

over electron Bessel functions which appear in the full dispersion relation. 

This implies that the validity of (2.1) is restricted to the domain 0 < 4° 

(with B =1.1) or e < 5° (for B =0.6). In contrast. Lemons and Gary [1977] 

retained the suras over Bessel functions and found that although Y(kjjj,6) was 

still maximum at small values of 9 , it remained relatively constant to 

beyond 9 < 20°. (Note that our definition of 0 is the complement of the 

definition used by Lemons and Gary [1977].) Figure 8 shows Yjj, and for the 

modified two-stream instability as a function of V./V. for B =1.1 and 0.6 

die 

(T = 10 and 5 ev, respectively). At threshold V.=3-^V. and tu << w, As 
e aim Ln 

V^ increases, so does Yj^, while 0 )^^^ increases more rapidly, eventually 

exceeding 
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As mentioned above, the scale of the MHD pulse can be estimated from 
(3.2), now rewritten in a more convenient form: 

At 1 AU we have fJ /m = 5.9x10~^. Assuming V./V. a 3 at threshold, c/V. = 

€ 0 u 1 0 

8.4x10’, so that Ax = 150(c/<j ) a 300 km for AB/B = 3 This value for Ax 

P© 

is consistent with observational and theoretical evidence of thicknesses of 
large amplitude MHD waves and shocks [see, e.g. Fredricks and Coleman, 1969 
and Manheiraer and Boris, 1972]. 

The situation at 2 AU is similar, and is also illustrated in Figure 8. 

Again the instability threshold is reached at V^/Vj^ = 3-4, and Ax = 92(c/aig) 

= 370 kra, where now n /m = 5x10"^, aB/B = 3, N=1.75 cra”^, (T.), = I .3 eV 

e e 1 

and Vgy = 450 kra s“^. This value for ax is close to the thicknesses of 

•^1000 kra reported by Smith and Wolfe [ 1977] for forward and reverse shock 

pairs using 453 kra s“^ as the average propagation speed. 

Saturation Effects . The modified two-stream instability results in 

heating both ions and electrons. The unraagnetized ions are primarily heated 

by raising while the electrons, constrained as thiey are to move along 

are primarily heated in (T ) , , [Ott et al., 1972], Stabilization of the 

MTSI has been studied using numerical simulations [Ott et al., 1972 and 

McBride et al. , 1972] where trapping was found to be the stabilization 

1 /2 

mechanism. Waves propagating at angles 6 > (m/M) =1.3° trapped elec- 

1 /2 

trons, while ion trapping occurred first if 6 < (m/M) . From Figure 7, 

significant growth is obtained only for 0 ^ < 1 . 3 °, so that ion trapping is 
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likely to be the more important effect. Consequently, (2.3) can be used to 
estimate the magnitude of the electric fields associated with this instabi- 
lity, and we find that for » 4 at both 1 and 2 AU, and E a 100 pV/m 

which increases to =300 .jV/ra when V./V, = 8-10. If the MTSI evolves to 

d i 

saturation, would be reduced to the marginally stable level (V^/V^ ® 3-4) 

and an amount of free energy equal to { 1 /2)Nm* ( AV^) ^ could be extracted from 

the current and converted to wave energy and heating [Lamps et al. , 1975]. 

The amount of electron heating that actually occurs depends sensitively on 

AV^; ranging from as much as a 10$ Increase in ® to a 

negligible change in (T„),, if aV. = V.. 

e 1 1 a 1 

The modified two-stream instability is not the only cross-field 
instability that can be excited in MHD pulses. Several others produce 
electrostatic waves at frequencies well above ar.u thus may be more 

easily observable by plasma-wave experiments. We will briefly consider two 
such instabilities; the beam-cyclotron and the ion-acoustic instabilities. 

Beam-cyclotron instability . Like the modified two-stream, the beam 
cyclotron instability (also known as the electron cyclotron drift instabi- 
lity) 1s driven by a relative drift of electrons and ions. In this case the 
instability arises from a coupling of Bernstein and ion acoutio waves 
[Forslund ct al. , 1970 and Gary and Sanderson, 1970]. The resulting 

unstable waves are confined to a narrow con? about the normal to £ (as are 
Bernstein waves) with frequencies close to harmonics of 0^, and kLg>1 [Gary, 
1971]. To excite the instability, V. must exceed C . Significant growth 
occurs even with T^ = T. and = 1 [Lashmore-Davies , 1971]. This instabi- 
lity has been extensively studied both analytically and in numercial 
simulations by Lampe et al. [1972]. 
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The dispersion relation for Js*£,=0 and kL >1 has been derived by Lampe 
et al. [1972], and is given by 

1 - /T, )-zY ^ + (k)i )“^ + (kX )"^ • 

® ^ \/2kV./ ® ® 2/2kV 

i e 

- 2/ ^ \ + i-cot(Wn^)-[z/„J:L_'\ + ?■( V t 

\/2kV^/ \/2kVg/ \/2kVg/ \/2kVe/ 

(3.4) 

where X^rv^/w^ is the electron Debye length, and is assumed parallel to 

Jl. The neglect of electromagnetic effects in (3.4) is well justified for 

the parameter range of interest [Lampe et al., 1972]. Because of the strong 

Landau damping of Bernstein waves propagating at 8^0, the unstable waves are 

similarly confined to angles e<n^/ (2irw k . X^) [Karaimura et al. , 1978]. At 1 

AU, (i)g/nga170, so that 0<1" for kx^=0.1. For a given frequency, unstable 

solutions to (3.4) are somewhat sensitive to the ratio of V./V. — large V, 

d i d 

results in smaller kX^, and hence a wider propag^^tlcn cone. 

To compare properties of this instability with the MTSI we solved (3.4) 
using solar wind parameters typical of 1 AU with B^=1.1. As in Figure 8, 
Vgw==450 km s , M=7 cm , B=5 ir, and Tg/T^-3. The results are shown in 
Figure 9 for V^/V^=10, where y and a) are plotted as functions of k for the 
first 8 harmonics of n^. The maxima in y, denoted y^, fall approximately at 

nn 

e 

- rw(i + k^x*)^''^ (3.5) 

a relationship derived by Lampe et al. [1972] for the case Tg>>T^. Although 
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the ratio is rather small, ranging from 5x10 - 5x10 , is 5-6 

times greater than found for the MTSI (of. Figure 8). On the other hand, 
the beam-cyclotron instability requires a slightly higher threshold for 
(5-6 V^, compared to 3-^1 for the MTSI). Thus it is not altogether clear 
which instability will dominate, although one can imagine a situation in 
which the MTSI is first excited, but is not immediately able to dissipate 
the available free energy. The MHD waves will then continue to steepen, and 
will increase until the beam-cyclotron instability is excited. Once 
excited, because of its large growth rate, the beam-cyclotron instability 
will saturate before the MTSI. The quasilinear and nonlinear stages of 
evolution have been examined by Lampe et al. [1972] who found that initially 
the beam-cyclotron instability saturated at a relatively low level via 
resonance-broadening. Once this happened, they found that the plasma could 
be further unstable to the ion-acoustic instability because the nonlinear 
dispersion relation for the beam-cyclotron instability when saturated 
through resonance- broadening has the form of the linear zero magnetic field 
ion-acoustic dispersion relation. In the solar wind, where T /T. is usually 

6 X 

too low to excite the ion-acoustic instability, the beam-cyclotron instabi- 
lity will saturate, allowing the modified two-stream instability to evolve 
until ion or electron trapping has reduced to a marginally stable state. 
However, because both the MTSI and beam cyclotron instabilities can heat 
electrons at least to some extent, it is quite possible that T will 
increase by just enough to excite the ion-acoustic instability. We explore 
this possibility below. 

In spite of the consensus that saturation of the beam-cyclotron 
instability via resonance broadening occurs at a low level of electric field 
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turbulence tf-ampe et al., 1972; Biskamp, 1973 and Lemons and Gary, 1976], it 

la rjot clear that this will be true for solar wind parameters . The ci cal 

amplitude of the turbulent fields when resonance broadening becomes impor- 

1 /? 

tant was found by Lamps et al. [1972] to be E ^ 96*[NT^(fl^/u^)Ml/kX^)] 

— ^ -2 

(raV/m), with T measured in eV, Assuming cm"^, kX “5 j< 10 , and w /O « 

170, we find E»1 mV/m, which is actually larger than the value of'X'100 pV/m 

found above for the saturation level of the modified two-stream instability. 

From our numerical solutions of (3»^), we know that the instability can 

be excited up to at least the eleventh harmonic of p , the exact maximum 

being a function of V^/V^ [Lampe et al., 1972]. These high harmonics have 

relatively large values of kX , and are thus confined to a very narrow cone 

0 

about 0=0 of much less than 1°. This may in turn limit the growth of the 
highest harmonics because of the necessity of the magnetic field to remain 
constant in direction over a growth time [Lemons and Gary, 1978]. These 
waves have frequencies in the range 100-1000 Hz, depending on the amount of 
Doppler shift, and should be easily observable. In fact, Wu and Fredricks 
[1972] have argued that this instability may have been observed in the 
earth's bow shock where amplitudes of 1-20 mV/m have been reported at 
frequencies near 1 kHz [Fredricks et al, , 1970; Fredricks and Coleman, 

1969 ]. Biskamp' a [1973] suggestion that the observed turbulence levels were 
too high to be explained by the beam-cyclotron Instability should be 
reexamined , at least for space plasmas . 

The quasilinear analysis and numerical simulations carrried out by 
Lampe et al. [1972] indicate that as long as V^>>Cg, the bulk of the energy 
taken out of the relative streaming between the prot;>na and electrons goes 
into heating electrons. Thus the beam-cyclotron instability appears to be 
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more efficient at heating electrons than is the modified two-stream instabi- 
lity, and with = ICV^ « 6C^, electron heating should be fairly efficient. 
At that drift velocity, only a small fraction of the total available free 
energy (equal to MV*) need go into electron heating in order that the 
electron temperature increase substantially. For example, would double 
if an amount of energy equal to MV^/30 were removed from the stream. Again, 
baaed on the results of Lamps et al. [1972], this is probably the maximum 
amount of energy that could be extracted before saturation by ion trapping 
occurs. Even this amount ■’3 unlikely unless the beam-cyclotron instability 
in turn excites the ion-acoustic instability. Therefore it is of interest 
to see Just how much it is necessary to heat the electrons before excitation 
of the ion-acoustic instability becomes possible. 

lQ.n--AGOustic Instability . The threshold for the ion-acoustic instabi- 
lity can be founu j i the dispersion relation discussed by Fried and Gould 
[1961] and Stringer [1964]. In our notation: 


1 - (1/2k*XM *(T^/T, )*2''f 

6 G X \ 


^m-kV 


/2kV, 


- (1/21;*X*) 
e 


= 0 


(3.6) 


For = 10, unstable roots first appear at T^/T^ » 10. Recall we have 
assumed that initially j j » 3* As we have seen, an increase in 
Tg/T^ by more than a factor of two is highly unlikely. If the ion-acoustic 
instability is not excited, relatively little streaming energy will be 
removed from the plasma. Therefore, the MHD pulse, or shock will continue 
to evolve, causing to increase still further. As approaches 20, it 
is then possible to excite the ion-acoustic instability with (T ),,/(T.),=6. 

® 1 1 i 1 
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The lon-acougtlc instability saturates by ion trapping [Biskamp and Chodura, 
1971 and Lamp® et al., 1972], and is capable of reducing the relative 
streaming to maginally stable levels. As an example, in Figure 10 we have 
plotted ui and y against k for the case V /V. = 20 and T /T. = 7. These 

Q X 6 X 

waves too will be subject to a Doppler shift in the solar wind, but should 

be observable by plasma wave instruments currently flown. [Note that 

because we have assumed that the protons are streaming through stationary 

electrons, the values of ui in Fig. 9 can exceed u. and u/k>C .] Saturated 

X s 

turbulent field intensities of order 1 mV/ra can be expected. Furthermore, 
the ion-acoustic instability produces waves propagating into a wide cone 
about 9=0. In fact, for V^/V^»20, 46=80°. 
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n. CONCLUSIONS 


In the p^eceeding sections we have discussed several kinetic instabili- 
ties which we argue are important in determining the microscale structure of 
the solar wind stream interaction region. The selected instabilities are 
not meant to comprise an exhaustive list. We have Ignored, for example, all 
of the lower-hybrid drift instabilities driven by strong gradients in 
density and magnetic field, primarily because their physics is similar in 
many respects to the magnetized ion-ion and modified two-stream instabili- 
ties [v., e.g. Lemons and Gary, 1978]. Another example of a related 
instability we did not discuss is the cross field, current driven ion- 
acoustic instability [Barrett et al. , 1972] because this shares many 
characteristics with the field free instability and MTSI . 

Perhaps the best way to summarize our discussion is to imagine follow- 
ing a stream interaction region as it evolves outward in the interplanetary 
medium. If one assumes that the fast stream originates in a coronal hole, 
while the slower material comes from the adjoining region, then close to the 
sun a shear interaction will predominate. Al.though we have not discussed 
any kinetic instabilities in this region, one can imagine several dynamical 
elfects which might be important. Por example, if the gradients in the 
shear become so large that the fast and slow ions are separated by only a 
few Larmor radii, then an electromagnetic ion cyclotron instability would be 
excited if the s.iear velocity exceeded ~2.7V^ [Eviatar and Wolf, 1968 and 
Goldstein and Eviatar, 1973]. On a larger scale the interface could become 
Kelvin-Helraholtz unstable. Eventually, however, the interface between the 
two streams will steepen as the faster flow begins to overtake the slower 
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(v. Fig. 1). At the interface the flow can be thought of as having two 
components, one along the interface, parallel to the magnetic field, the 
other across the field. Because of their large gyroradii, the ions will be 
the first component of the fast stream to penetrate the slow plasma. 
Depending primarily on the relative velocities of the two streams and the 
velocity of the solar wind i.^. the slow flow, this situation will become 
unstable to the magnetized ion-ion instability somewhere between 1 and 2 AU. 
The magnetized ion-ion instability will heat the ions (proportional to their 
mass) . The amplified electric fields produced by the interpenetration of 
the two streams can then give ri^o to drifts which excite either the 
modified two-stream instability or cne of the other instabilities discussed 
in §2. Excitation of these kinetic instabilities will quickly erode any 
sharp gradients at the interface, leaving it in a state of marginal 
stability. 

Beyond one AU the compression region , which has gradually steepened as 
the streams evolve, forms pulses within which strong gradients in £. can 
excite various electron-ion cross field instabilities. The possible range 
of frequencies excited extends from the lower-hybrid (due to the modified 
two-stre&m instability) through the ion-acoustic frequency range, to 
harraon^ ' of the electron cyclotron frequency. It appears unlikely that 
either the modified two-stream instability or the beam-cyclotron instability 
would be able to reduce the streaming to a marginally stable state. Should 
that be the case, the compression regions would continue to evolve until 
those two Instabilities had heated the electrons sufficiently to drive ion- 
acoustic waves unstable. Once excited, the observed shock thickness and 
electron temperature would remain close to their marginally stable level 
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[Manheimer and Boris, 1972]. As yet, insufficient research has been 
conducted on the detailed structure of forward and reverse shocks in stream 
interaction regions to know whether the structure is consistent with this 
scenario. Hf'vever, Morse and Greenstadt [1976], in an analysis of the 
earth's bow shock, did conclude that under certain conditions the thickness 
of the shock appear to be controlled by the conditions for marginal stabili- 
ty against the ion-acoustic instability. Present and future plasma wave 

experiments on deep space prober should be able to answer many of these 
questions . 
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APPENDIX 


The dispersion relation can be written in the form: 


=0 (A1) 

where £ is the dyadic given by Montgomery and Tidman [1961 eq. 10.35]. It 
is convenient to write the electric field as the sura of electrostatic and 
electromagnetic waves: 

K = -itflijil + (m/c) A] (A2) 


where « and A ai’e the scalar electrostatic and vector electromagnetic vector 
potentials, respectively. We take in the 1 direction, Jc = (kj^,0,kj|), 
and A = Ax. Expansion of (A1) gives 




+ « A (o)/c) k, . (Dy^ •*- D^y) = 0 


(A3) 


where we have used D = -D . 

xy yx 

Following McBride et al. [1972], we assume fl,<<jtu|<<|n I, k,,«k,, 

i e ‘ I [ _L’ 

k*Lg^<<1, and k^L^^>>l, which implies that the ion trajectories are well 
approximated by straight lines and the electron trajectories by helicies. 
We also assume that the temperatures of each species are the same in both 
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streams and are isotropic. However, the ion and electron temperatures may 
differ from on another. 

For the magnetized ion-ion instability we take the following distribu- 
tion functions 

f, = J a exp[-(v - ii)*/2Vh + (1-a) expL-v*/2Vn 

(2n)3/2v3 

f’e = ^ exp(-v^/2VM (A4) 

(2v)3/2v. 

e 

where ii = (Usinr ,0 ,Ucosr) is the drift velocity and r is the garden-hose 
angle . 

Let Ip be the angle between Js and ii, and 0 the angle between Jt and the 
normal to Then the ion component of the electrostatic dispersion 

relation, which is the part of (A3) multiplying 4>*, becomes 

1 - ^ a I' h " + (1 - a) z V 

2k^V| \ /2kV^ ) 

where Z(^) = (1 //ti) dy exp(-y^)/(y - c) is the plasma dispersion 

function. Using ip + 0 + r = ir/2, and sin0 = 0, we have 

Jc*Ii = k U Gosip = k U (sinr cos0 + cosT sinO) 

= k U sin( r + 0 ) = k U sinr cose. 

This provides the ion terms of equation (2.1). 
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The electrons, on the other hand, are highly magnetized. In the strong 
magnetic field limit, the electrostatic electron terms become 

m * 

__L {1 - exp(-k*LM I„(k'L *) [1 + (1/2)Z'(u//2kV sinQ)]) 

e o e e 

e 

Iq is the modified Bessel function. 

The electroraangetic and mixed terms proportional to A* and $A in (A3) 
tend to stabilize these instabilities, and so must be kept [McBride et al. , 
19723. However, in the range of frequencies of interest here, it is 
possible to further simplify these terras by taking the fluid limit in which 
all resonances are ignored. Thus we take the arguments of the Bessel 
functions to be small, the arguments of the plasma dispersion functions to 
be large (except for n = 0) and expand accordingly. We also use 1^ = Ij, 

= 1 + (1/^) *(kV^/fi^) and I, = (1/2)* (kV /n )^. We retain the warm 

electron terms for reasons discussed in the text and obtain terms of the 
type 


(to /kc) * 

1 + f 

1 + + w /kc 

e e 


Z^(u)//2kV sine) + (m /a )* 
e e e 


1 + 


tD^/kc 
1 + B 


e 


Combining these terms gives (2.1). If we replace the ion distribution 
function given in (A4) by 


f 


i 


1 

(2,r)3/2v» 


exp[-(i - Ijj) V2V?] 


we obtain the dispersion relation for the modified two-stream instability 
given in (3.1). 
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FIGURE CAPTIONS 


Fig. 1 . Schematic view of the stream interaction region showing the compres- 
sion (shaded) and rarefaction produced as high velocity plasma (long arrows) 
overtakes the low velocity plasma. The view is onto the ecliptic plane from 
above the north pole. The spiral pattern of the magnetic field is also 
shown. For further details see Pizzo [1979a], from which this figure is 
adapted. 

fig. 2. Growth rate y, and real frequency u, of the •magnetized ion-ion 

instability (eq. 1.1) at 1.5 AU, In panel (a), tu and y are plotted against 

0, the angle between and the normal to £, v;ith |iij held constant, and 

equal to the value which maximizes y(k). In panel (b), the wavenumber 

variation of w and y are shown with 6=6^^, the value which maximizes y(6). 

'K'te that the maximum value of y is relatively insensitive to variations in 

B , although the angular spread of y i.T greatly reduced as B decreases. 
6 e 

Because w is relatively insensitive to variations in B , only iii(B =0.9) is 

e e 

plotted. In computing these values of oj and y we used Vg^= 300km s~^ and 1 
AU values of T^^ and N equal to 2 eV and 7 cm“^, respectively. The scaling 
laws defined in the text were employed to find T^ and N at 1 .5 AU. In 
addition, we used U=40 km s \ where Usinr is the component of the velocity 
jump in Vgy at the interface that is perpendicular to 

Fig. 3. The dashed curve represents a solution to (1.1) for the magnetized 
ion-ion instability in wliich the electron plasma dispersion function was 
approximated using the asymptotic expansion valid for large argument, or 
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equivalently, small electron thermal velocity. That approximation is seen 

to be valid only for very small angles 3t least for parameters 

typical of the solar wind. In this example, the 1 AU values of U, Vg^, Tj|^ 

“1 -1 

and 0 were 50 km s , 300 km s , 1.5 eV and 1.9| respectively. The 

solution was obtained at R=2Alf. T.(1 AU) and 0 (1 AU) were found from the 

1 e 

assumed radial dependence discussed in the text. 

Fig. A plot of ID and as functions of heliocentric distance, with 

mm ’ 

* 1 

Vgj^=300 km s . Yjjj is the maximum growth rate of the magnetized ion-ion 
instability, maximized with respect to both k and 6; is the accompanying 

value of (j(k ,0 ). Results are plotted using both U=40 and 50 km s”^. Also 
plotted is Uq(R), given by (1.2), and n^(R). Note that Justifying 

the assumption that the ions are unmagnetized. The largest values of Yjjj are 
found between 1 and 2 AU. 

Fig. 5 . The threshold conditions of the magnetized ion-ion instability 

plotted as a function of U and R for three different values of Vg^. The 

magnetized ion-ion instability is unstable for values of U which lie above 
the curve corresponding to the choice of 

Fig. 6. The maximum value of the electric field in lower-hybrid waves 

excited by the magnetized ion-ion instability. Stabilization is assumed to 

occur by ion trapping. The largest values of E are found between 1.5-2. 5 

-1 -1 

AU. The curve was obtained using Vgy=300 km s , and U=50 km s 

Fig. 7 . In panel (a), Y(k,0^), the growth rate of the modified two-stream 


AO 



inotability , is plotted as a function of k for Bg=0.£ and 1.1 at 8 = 0jj,* In 
panel (b) , k is held constant at k=k^, and y is plotted against e. In both 
panels Vg^=iJ50 kra s \ R= 1 AU, T^=3-3 eV, 3=5 y, and N=7 cm”^. In addi- 
tion, to(k,0 ) and u(k ,0) are shown for 6 =1.1. 
ram e 


Fig. 8. Real and imaginary frequencies of the modified two-stream instabili- 
ty as a function of V^/V^. Solid lines represent the solution at 1 AU using 
N=7 cm and (T^)j^=3.3 eV. Dashed lines represent the solution at 2 AU 
using N=1.75 cm"^, and (T^)^=1.3 eV. In both oases Vgj^=H50 km s"^. The 

solution for a>„ is shown for B =1.1 at 1 AU and 6 =1.5 at 2 AU, while y is 

m e e ’ 'm 

also plotted for 6^(1 AU)=0.6 and 8^(2 AU)=O.0. The variation in N, T, , and 

6 6 1 

8g between 1 and 2 AU is a result of the scaling with heliocentric distance 
discuessed in §1. 


Fig. 9 . Real and iraa*lnary roots of the dispersion relation for the beam- 

cyclotron instability as functions of kL. (bottom scale), and kA„ (top 

6 6 

scale). The parameters are characteristic of 1 AU (N=7 om“^, B=5 y, 

T /T, = 3, B =1.1). In addition, V^„=i<50 km s"^ and V^/V, = 10. 
e 1 e &W d i 

Fig. 10. Solution of the ion-acoustic dispersion relation assuming V^/V^=20 
and T /T.=7. As before N=7 cm”^ and V„,,=l150 km s~\ 
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Figure 3 
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